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NATO STO Task Group AVT-201 on “Extended Assessment of Reliable Stability & 
Control Prediction Methods for NATO Air Vehicles” is studying various computational 
approaches to predict stability and control parameters for aircraft undergoing non-linear 
flight conditions. This paper contributes an assessment through correlations with wind 
tunnel data for the state of aerodynamic predictive capability of time-accurate RANS 
methodology on the group’s focus configuration, a 53° swept and twisted lambda wing 
UCAV, undergoing a variety of roll, pitch, and yaw motions. The vehicle aerodynamics is 
dominated by the complex non-linear physics of round leading-edge vortex flow separation. 
Correlations with experimental data are made for static longitudinal/lateral sweeps, and at 
varying frequencies of prescribed roll/pitch/yaw sinusoidal motion for the vehicle operating 
with and without control surfaces. The data and the derived understanding should prove 
useful to the AVT-201 team and other researchers who are developing techniques for 
augmenting flight simulation models from low-speed CFD predictions of aircraft traversing 
non-linear regions of a flight envelope. 
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wing span, =1.5 3 8m 
reference span (b/ 2), =0.769m 
lift coefficient, =Lift/q CXj S re f 

rolling moment coefficient about x-body axis, =Rolhng_Moment/goo5W/£re/ 

pitching moment coefficient about y-body axis, =Pitching_Moment/^ 0 o5 , re /c re / 

yawing moment coefficient about z-body axis, =Yawing_Moment/g 00 Sre/&re/ 

normal force coefficient, =Normal_Force/#o 0 ‘SW/ 

pressure coefficient 

wing root chord, = 1.061m 

reference chord, =0.479m 

sinusoidal oscillation frequency about pitch, roll, or yaw axis, Hz. 
reduced frequency for pitch oscillation, =2jzf- c re f/Uao 
reduced frequency for roll and yaw oscillation, =2jtf • (b/2) /Uao 
left-outboard/left-inboard trailing-edge control surface deflections, positive downward 
order of magnitude drop in solution residual, e.g. -4 means 4 orders of magnitude drop 
freestream Mach number 

Moment Reference Point for pitching moment, =0.600m from apex 
Number of time steps per oscillation cycle, =2jt/(k 0 At ) 

Dynamic pressure, N/m 2 
Reynolds number based on c re f 

right-outboard/right-inboard trailing-edge control surface deflections, positive downward 

DNW-NWB rear rotation point for pitch and yaw oscillation, =0.855m from apex 

primary& secondary VGRID “viscous” stretching factors, see Eq. 1 

reference area, =0.7 7m 2 

physical time step, seconds 

characteristic time step, =At 0 UJc re f 
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= freestream velocity, ~50m/s 
= angle of attack, deg. 

= angle of sideslip, deg. Positive nose left (wind from right). 

= VGRID “viscous” grid spacing normal to surface at node i+1, see Eq. 1, meter 
= spacing of first node off of surface in “viscous” grid layers, see Eq. 1, meter 
= nominal roll angle, deg. 

= range of roll oscillation about <P 0 through body axis, +/-deg. 

= nominal pitch angle, deg. 

= range of pitch oscillation about <9^ through body axis, +/-deg. 

= nominal yaw angle, deg. 

= range of yaw oscillation about *7^ through body axis, +/-deg. 
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Applied Vehicle Technology 

Computational Fluid Dynamics 

DFR-F 1 9 with undeflected control surfaces 

DFR-F19 with deflected control surfaces, FOB/FIB=-20/-20, ROB/RIB=+20/+20 deg 
Degree-of-freedom 

German-Dutch Wind Tunnel located in Braunschweig 
Forced-oscillation sinusoidal motion 


KE, k-e 

NATO/STO 

RANS 

SA 

SID 

SST 

UCAV 


Jones-Faunder linear k-epsilon two equation turbulence model 

North Atlantic Treaty Organization / Science and Technology Organization 

Reynolds Averaged Navier Stokes 

Spalart-Allmaras one equation turbulence model 

System identification 

Menter’s Shear Stress Transport two equation turbulence model 
Unmanned Combat Air Vehicle 


I. Introduction 

Historically, many civil and military aircraft development programs have encountered costly stability and control 
(S&C) deficiencies during early stages of flight test despite thousands of hours of wind tunnel testing. These 
surprises have occurred across the speed range from takeoff and landing to cruise flight, and particularly at the 
fringes of the flight envelope where separated flows dominate. It is widely believed within the research community 
that the next significant advancement in the state of the art for predicting the S&C characteristics of a new vehicle 
might be through the application of computational fluid dynamics (CFD) tools. The application of high-end CFD 
codes with Reynolds-averaged Navier Stokes (RANS) or better level of technology to specific areas of S&C interest 
before first flight can help focus wind tunnel programs and provide improved understanding of the underlying flow 
physics. A focus for such an advancement has been fostered by a coordinated international team of researchers 
under the NATO/STO AVT-161 Task Group titled ’’Assessment of Stability and Control Prediction Methods for 
NATO Air & Sea Vehicles” (2007-2010) [1] and the AVT-201 Task Group titled “Extended Assessment of Reliable 
Stability & Control Prediction Methods for NATO Air Vehicles” (2011-2015) [2]. Under AVT-201, the focus has 
been extended to include investigation of applying system identification (SID) techniques to CFD-derived S&C data 
for constructing reduced-order models for flight simulators. NASA is also leveraging this NATO activity to guide a 
parallel study of S&C CFD modeling and SID for civil transports near stall within the Vehicle Systems Safety 
Technology project under the NASA Aviation Safety Program [3]. 

The prevalent experimental methodology for quantifying the dynamic stability of an aircraft calls for extracting 
the dynamic stability derivatives from the force and moment coefficients as the aircraft undergoes a periodic motion 
about a body axis in a wind tunnel. Experimentalists within the AVT-161 and AVT-201 teams have conducted 
dynamic wind-tunnel tests on the focus configuration, called Stability and Control Configuration (SACCON), 
undergoing pitch, roll, and yaw oscillation to measure the dynamic stability and aerodynamic characteristics [4,5]. 
Other team members are running parallel dynamic computational studies following a similar approach of modeling 
SACCON undergoing pitch, roll, and yaw oscillation with their respective CFD tools [6-12]. Alternative approaches 
applying system identification techniques to multiple frequency/amplitude motions [13], and more recently 
exploring indicial maneuvers [14] that have the potential of reducing the computational cost for constructing 
dynamic derivatives by a factor of five. The collective experiences from these studies have confirmed the difficulty 
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of achieving accurate and efficient computation of dynamic stability characteristics for the highly nonlinear portions 
of the flight envelope. 

For purposes of coordinating a focused international research study, the AVT-201 Task Group defined a set of 
common static and dynamic SACCON test cases for team members to compute [15]. These cases include static 
angle of attack (alpha) and sideslip (beta) sweeps and select dynamic forced-oscillations about the pitch and yaw 
axes for various combinations of control surface deflections. This paper documents a NASA CFD investigation of 
SACCON S&C characteristics for subset of the AVT-201 common test cases. The roll oscillation solutions are 
computed for an AVT-161 test case without control surface deflections. A limited static study is performed to gain 
understanding of grid and turbulence model sensitivities, which will provided guidance for dynamic force- 
oscillation computations that follow. 


II. Geometry and Experiment 

A. General description and reference conditions 

The SACCON geometry presented in Fig. 1 is a generic representation of a UCAV configuration with 5 3 -deg 
swept leading-edge. The variation used in this investigation has a sharp inboard leading edge segment that 
transitions to a “medium” round leading edge on the outer wing panel. This outer panel has parallel leading and 
trailing edges that twist 5-deg in washout. As noted in Fig. 1, the wingspan is 1.54m, root chord is approximately 
lm, and reference chord, c re f is 0.48m. The reference area, S re f is 0.77m 2 . For pitch and yaw forced oscillation 
testing, the rear rotation point ( RRP ) is positioned near the aft end of the root chord, downstream of the aerodynamic 
moment reference point, MRP. With roll oscillation testing, the model rotates around an aft-mounted sting through 
the longitudinal axis. 

Two SACCON wind tunnel models were used to provide correlation data for this study; the DLR-F17 used as 
the primary test article in the prior AVT-161 task group, and a new model with control surfaces designated DLR- 
F19. Both were constructed with the dimensions noted in Fig. 1 from lightweight composite material to facilitate 
dynamic oscillation testing. Both were equipped with pressure taps on the upper and lower surfaces of the model 
that are configured for dynamic measurement of unsteady pressures. 

Initial tests with the DLR-F17 SACCON configuration during AVT-161 revealed an arbitrary boundary layer 
transition line on the upper surface of the model detected by infrared thermography. This discovery led to a decision 
to apply carborundum trip grit to fix transition near the leading edge. Transition of the boundary layer to fully 
turbulent flow was subsequently verified with this configuration mounted on a belly sting [4]. 

B. DLR-F19 SACCON with belly sting 

The DLR-F19 wind-tunnel model was mounted to a six-component balance by a post sting attached to the underside 
at the RRP. The model is shown in Fig. 2 (left) mounted inverted in the closed test section of the German-Dutch 
Wind Tunnel (DNW-NWB) located in Braunschweig, Germany with the control surfaces deflected in the 
LOB/LIB=-20/-20, ROB/RIB=+20/+20 deg configuration. The location of the pressure rows on the DLR-F19 is 
shown in Fig. 1. A yaw-link support sting with 15° crank angle was used for the static data and yaw oscillation 
cases. The pitch oscillation runs utilized a pitch-link support sting with 6° crank angle and a 7th axis push rod. 
Forced motion was induced using the DNW-NWB Model Positioning Mechanism (MPM) that is a six degree-of- 
freedom (DOF) parallel kinematics system designed for static as well as for dynamic model support. The MPM is 
located above the test section and can be operated in the open test section as well as in the closed one. The location 
of oscillation axes can be chosen arbitrarily and in addition to classic sinusoidal oscillations the MPM can perform 
multi-DOF maneuvers. Details of the model, MPM, and the installation are described in Ref. [4]. 

In future references to the DLR-F19 configuration in this paper, the cases with undeflected control surfaces will 
be designated CSO. The cases with deflected control surfaces LOB/LIB=-20/-20, ROB/RIB=+20/+20 deg will be 
referenced as CS20. 

C. DLR-F17 SACCON with aft sting 

The DLR-F17 wind-tunnel model is shown in Fig. 2 (right) supported by an aft sting mount in the NASA 
Langley Research Center (LaRC) 14- by 22-Foot (4.3- by 6.7-meter) Subsonic Tunnel in Hampton, Virginia, United 
States. Forced sinusoidal motion was applied about the roll axis from a motor on the sting support system. The 
DLR-F17 was the primary SACCON test article in the AVT-161 task group. Hence, extensive force and moment, 
pressure, and Particle Image Velocimetry (PIV) flow measurement static and dynamic data is available for this 
configuration mounted on both belly sting and aft sting [4,5]. The surface of the model was painted with shiny black 
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paint that contains particles of Rhodamine B to reflect laser light at a different wavelength, which can be filtered so 
as to produce highly accurate PIV measurements close to the surface. 

III. Methodology 

The CFD computations are performed using the NASA Tetrahedral Unstructured Software System (TetrUSS) 
[16]. This system consists of loosely integrated, user-friendly software that comprises a geometry setup utility 
GridTool, a tetrahedral grid generator VGRID, a flow solver USM3D, and post-processing visualization and data 
extraction utilities. Body motion is controlled through a third-party process control framework called FD-CADRE. 

A. Grid generator 

VGRID [17,18] is a tetrahedral grid generator based on the advancing front method (AFM) for generation of 
surface triangles and “inviscid” field cells, and the advancing layers method (ALM) for generation of thin-layered 
“viscous” cells. Both techniques are based on marching processes in which tetrahedral cells grow from an initial 
front (triangular surface mesh) until the volume around the geometry is filled. Unlike the conventional AFM, which 
introduces cells into the field in a totally unstructured manner, the ALM generates organized layers of thin 
tetrahedral cells, one layer at a time, while maintaining the flexibility of AFM. Once the advancing front process is 
completed in VGRID, an additional postprocessing step is required using POSTGRID to close any open pockets and 
to improve grid quality. VGRID input files are generated by an interactive geometry manipulation program, 
GridTool. This graphics tool can import surface definitions from initial graphics exchange specification (IGES) files 
containing nonuniform rational B-spline surfaces and curves, as well as PLOT3D point definition files. GridTool is 
used to manipulate the geometry and to define necessary geometric surface patches and grid spacing (source) 
parameters. 

B. Flow solver 

The computations are performed with USM3D [19], which is a parallelized tetrahedral cell-centered, finite 
volume Navier- Stokes flow solver. The term “cell centered” means that the finite volume flow solution is solved at 
the centroid of each tetrahedral cell. Inviscid flux quantities are computed across each tetrahedral cell face using 
various upwind schemes. Spatial discretization is accomplished by a novel reconstruction process, based on an 
analytical formulation for computing solution gradients within tetrahedral cells. The solution can be advanced in 
time by a 2 nd -order “physical” time step scheme, a 2 nd -order “dual” time step scheme, or to a steady-state condition 
by an implicit backward-Euler scheme. Several turbulence models are available: the Spalart-Allmaras (SA) one- 
equation model, the two-equation k - 8 turbulence model, the Menter Shear Stress Transport (SST) two-equation 
model, and the nonlinear Algebraic Reynolds Stress Models of Girimaji and Shih/Zhu/Lumley. Detached Eddy 
Simulation (DES) has been implemented in all of the turbulence models, but has not been fully tested. The latest 
extensions to the USM3D flow solver are described in Ref. [20]. 

C. General motion controller 

The USM3D flow solver has been installed as a plug-in to the FD-CADRE framework (Fluid Dynamics - 
Computational Analysis of Dynamically Responsive Environments) [21] developed at Arnold Engineering 
Development Center in Tullahoma, Tennessee. FD-CADRE is a generalized dynamic process control manager for 
coupling various plugins, e.g. flow solver, 6-DOF motion generator, aeroelastic structural module, etc. The 
advantage is that one can couple an existing legacy code such as USM3D to additional tools that dramatically 
expand its analysis capabilities. The USM3D/FD-CADRE combination has been used heavily by users in the 
Department of Defense for computing complex moving-body problems with overset tetrahedral grids. The current 
work only utilizes the prescribed motion capability on non-deforming, non-overset grids. 

D. Solution strategy 

USM3D solutions are generated using an implicit second-order physical time-step scheme for both the static and 
dynamic cases. Inviscid fluxes are computed with Roe’s flux difference scheme without limiting. The time-accurate 
RANS flow solutions are primarily computed with the SA turbulence model. Limited assessments are also made for 
the static cases with the SST and k-8 turbulence models. The boundary conditions consist of a surface no-slip 
constraint on the wing and support sting and characteristic inflow/outflow on the outer box. The forced oscillation 
(F-O) solutions are generated using nondeforming solid-body rotation of the full grid about the x, y, or z axis and 
initialized by restarting from a converged static solution at the prescribed angle of attack, a 0 , and cycled for two full 
sinusoidal oscillations. However, the hysteresis of forces and moments is converged to its periodic solution after the 
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first quarter cycle of oscillation; only the second cycle is plotted in the correlation with experiment. Details of 
specific time step and turbulence model strategies will be included in Section VI. 

IV. Grid Generation 

A series of full-span tetrahedral grids were generated for the DLR-F 1 9/F 1 7 S ACCON configurations shown in 
Fig. 3 using a developmental version of VGRID [18] and following guidelines developed in Ref. [6]. Thin-layer 
tetrahedral grids were generated to meet requirements for cell-centered computations from the USM3D flow solver. 
A near-wall first-cell spacing was prescribed, based on flat-plate turbulent boundary layer theory, to achieve a 
tetrahedral cell centroid turbulent wall coordinate (y + ) of 0.5 at a longitudinal distance of 0.5*c re /for a Re cre f= 1.5 
million. Since the VGRID Advancing Layers Method marches nodes away from the vertices of the surface triangles 
(which are subsequently connected to form tetrahedral cells), an initial VGRID spacing, Si, corresponding to a y + =2 
at the first node was prescribed in order to achieve the y + =0.5 at the first cell centroid. Subsequent USM3D 
computations confirmed that an average first-cell y + ~0.4 was achieved. For the S ACCON at the nominal wind 
tunnel chord Reynolds number of 1.5 million, the required VGRID first-node spacing is d 1 /c re j=2.4E-05, and 
stretching factors of ratel=0.15 and rate2=0.02 , where the nodal spacing layers are defined by the Eq. (1). 

S i+1 = -[1 + rate 1 ■ (1 + rate 2)7 (1) 

The same “viscous” spacing distribution was applied to all DLR-F 19 and DLR-F 17 grids and resulted in 
approximately 72 tetrahedral cells (24 nodes) across the boundary layer at the mid-chord of c re f. 

Six grids were generated on the DLR-F 19 configuration with the 15° yaw-link post sting; coarse/medium/fme 
(G1/G2/G3) with and without flap deflections with the characteristics listed in Table 1. Similarly, 
coarse/medium/fme grids were generated on the DLR-F 17 with aft roll sting producing the characteristics listed in 
Table 2. In both cases, the VGRID sources from the grids in Ref. [6] were transferred into the VGRID input files for 
the current configurations, thereby insuring a consistency with the prior study. 


Table 1 - General full-span mesh properties for DLR-F19 with 15° yaw-link belly sting. 


Grid 

LOB/LIB 

[°] 

ROB/RIB 

[°] 

VGRID 

ifact 

Tets 

[millions] 

Nodes 

[millions] 

Wing 

triangles 

Wing 

nodes 

G1.00 

- 

- 

1.9440 

7.233 

1.228 

73,456 

36,738 

G2.00 

- 

- 

1.4215 

12.967 

2.203 

132,474 

66,251 

G3.00 

- 

- 

1.0000 

25.776 

4.383 

261,050 

130,543 

G1.20 

-20/-20 

+20/+20 

1.9440 

10.286 

1.749 

120,364 

60,192 

G2.20 

-20/-20 

+20/+20 

1.4215 

17.963 

3.056 

212,456 

106,242 

G3.20 

-20/-20 

+20/+20 

1.0000 

34.960 

5.952 

414,984 

207,510 


Table 2 - General full-span mesh properties for DLR-F17 with aft-roll sting. 


Grid 

VGRID 

ifact 

Tets 

[millions] 

Nodes 

[millions] 

Wing 

triangles 

Wing 

nodes 

Gla 

1.9440 

6.523 

1.228 

71,706 

35,881 

G2a 

1.4215 

11.689 

1.987 

130,048 

65,059 

G3a 

1.0000 

23.185 

3.943 

256,664 

128,379 


The upper surface grid triangulations for G1.20, G2.20, and G3.20 are shown in Fig. 4 for the DLR-F 19 with 
deflected control surfaces (CS20), along with corresponding volume grids slices through the right-wing control 
surface in Fig. 5. The grid refinement was achieved by changing the “if act” parameter in the VGRID input, which 
globally scales the spacing for the “inviscid” portion of the grid, as highlighted in Fig. 5. The normal distributions 
across the “viscous” layers were not scaled in the process. Spanwise grid stretching with a ratio as high as 10-to-l 
was applied along the leading- and trailing-edges for each Gl, G2, and G3 as described in Ref. [6] in order to reduce 
cell count. Additional detail of the surface and volume grids in vicinity of the outboard gap for the right-wing 
control surface is provided in Figs. 6 and 7. The gap between abutting flap/wing surfaces was sealed. 
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V. Scope of Investigation 

For the AVT-201 effort described in the Introduction, the task group developed a set of common test cases for 
the DLR-F19 SACCON configuration to be analyzed and compared amongst the team. These included both static 
alpha and beta sweeps, and some cases for dynamic forced oscillation about pitch and yaw axis, all with various 
combinations of control surface deflection [15]. A subset of the static a- and /3 - sweeps will be used to assess 
sensitivities of the flow solution to grid density and turbulence models. The static results then provided guidance for 
pursuing the dynamic cases. The roll oscillation study is performed on the DLR-F17 SACCON for correlation with 
NASA Langley Research Center 14X22 subsonic wind tunnel data obtained by the AVT-161 task group. 

A. Static a and sweeps 

As depicted in Table 3, the static investigation includes two angle-of-attack sweeps from a=0-30 deg, one with 
and one without control surface deflection. Solutions are provided at cc= 0.0, 5.3, 10.0, 12.5, 14.7, 15.7, 16.8, 17.8, 
18.9, 19.9, 22.5, 25.0, 30.0 deg and zero sideslip. Static sideslip sweeps of / 3=+/-10 deg at A/3=l deg increments are 
performed at a= 10.48 and 14.67 deg. The corresponding correlation dataset from the DLR TN2545 wind tunnel test 
[22] are listed in the first column. All cases will be computed on the Gl, G2, and G3 grids from Table 1. The 
Spalart-Allmaras (SA), Mentor Sheer Stress Transport (SST), and k-8 (KE) turbulence models will be assessed on 
select grids. 


Table 3 - Static a/g-sweeps. SACCON DLR-F19. 


DLR 

TN2545 

LOB/LIB 

[°] 

ROB/RIB 

[°] 

a 

[°] 

p 

n 


R^cref 

[million] 

1001 

- 

- 

0-30 

0 

0.146 

1.57 

1103 

-20/-20 

+20/+20 

0-30 

0 

0.146 

1.57 

1007 

_ 

_ 

10.48 

+10...-10 

0.146 

1.57 

1109 

-20/-20 

+20/+20 

10.48 

+10...-10 

0.134 

1.46 

1008 

- 

- 

14.67 

+10...-10 

0.146 

1.57 

1110 

-20/-20 

+20/+20 

14.67 

+10...-10 

0.134 

1.46 


B. DLR-F19: Forced oscillation in pitch 

The AVT-201 common test cases for forced oscillation in PITCH are listed in Table 4. Two reduced frequencies 
for sinusoidal pitch oscillations (k P =0.06 and 0.12 ) are assessed for AO=+/-4.7 deg about the nominal pitch angle 
®0=1 0.3 deg, with and without control surface deflections. The corresponding correlation dataset from the DLR 
TN2545 wind tunnel test [23] are listed in the first column. The pitch oscillations were produced by a pitch link 
support with 6° crank angle and an axis push rod. Hence, there was relative motion between the SACCON and the 
post support sting that was not modeled in the computations. The Mach and Reynolds numbers are adjusted in the 
computations to correspond with their respective test conditions. The appropriate grids were selected from the 
results of the static study. 


Table 4 - Common dynamic test cases for PITCH F-O. SACCON DLR-F19. 

0 o =O°, W 0 =0°. 


DLR 

TN2545 

LOB/LIB 

[°] 

ROB/RIB 

n 

0o 

n 

AG 

[°] 

/ 

[Hz] 

kp 

2jtfc re JU„ 


Re cref 

[million] 

2342-50 

- 

- 

10.3 

4.7 

1 

0.06 

0.146 

1.57 

2351-59 

- 

- 

10.3 

4.7 

2 

0.12 

0.146 

1.57 

2666-74 

-20/-20 

+20/+20 

10.3 

4.7 

0.94 

0.06 

0.138 

1.49 

2675-83 

-20/-20 

+20/+20 

10.3 

4.7 

1.88 

0.12 

0.138 

1.49 


C. DLR-F19: Forced oscillation in yaw 

The AVT-201 common test cases for forced oscillation in YAW are listed in Table 5. Two reduced 
frequencies for yaw oscillation ( k Y =0.10 and 0.19) are assessed for AW=+I - 5 deg sinusoidal oscillations at a pitch 
angle of 6>o=10.1 deg, with and without control surface deflections. The run numbers from the corresponding 
correlation DLR TN2545 wind tunnel [23] dataset are listed in the first column. The Mach and Reynolds numbers 
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are adjusted in the computations to correspond to their respective test conditions. The appropriate grids were 
selected from the results of the static study. 


Table 5 - Common dynamic test cases for YAW F-O. SACCON DLR-F19. 

<Po=0°, W 0 =0°. 


DLR 

TN2545 

LOB/LIB 

[°] 

ROB/RIB 

n 

00 

[°] 

AW 

[°] 

/ 

[Hz] 

ky 

Moo 

R dcref 

[million] 

2270-78 

- 

- 

10.1 

5 

1 

0.10 

0.146 

1.57 

2279-87 

- 

- 

10.1 

5 

2 

0.19 

0.146 

1.57 

2162-70 

-20/-20 

+20/+20 

10.1 

5 

0.94 

0.10 

0.138 

1.49 

2171-79 

-20/-20 

+20/+20 

10.1 

5 

1.88 

0.19 

0.138 

1.49 


D. DLR-F17: Forced oscillation in roll 

The ROLL oscillation cases for the AVT-161 DLR-F17 configuration are listed in Table 6. These consist of 
prescribing sinusoidal roll oscillations of A<P=±5 deg at 7 frequencies for O 0 =10 and 20 deg. The corresponding run 
numbers from the NASA dataset T134 [4] are listed in the first column. The T134 Mach and Reynolds numbers are 
M^=0.053 and Re cre /=0.375 million. 


Table 6- Dynamic test cases for ROLL F-O. 
SACCON DLR-F17. & 0 =0 o , W 0 =0°. 


NASA 

T134 

O 

§ L- 

AW 

[°] 

/ 

[Hz] 

kR 

2n fbre/Uoo 

17 

10, 20 

5 

0.24 

0.063 

18 

10, 20 

5 

0.36 

0.095 

19 

10, 20 

5 

0.44 

0.116 

20 

10, 20 

5 

0.55 

0.145 

22 

10, 20 

5 

0.66 

0.174 

25 

10, 20 

5 

0.86 

0.227 

12 

10, 20 

5 

1.00 

0.264 


VI. Results and Discussion 

Correlations of computational and experimental results are presented to assess CFD predictive capabilities for 
the static and dynamic S&C characteristics of the SACCON configuration with and without control surface 
deflections. A cursory a-sweep is generated with the DLR-F19 to provide overlap with prior SACCON 
computational results in Ref. [6] and to assess longitudinal aerodynamics of the configuration with deflected control 
surfaces. Then a series of /Lsweeps are performed to assess grid and turbulence model sensitivities, and to 
investigate static lateral/directional S&C characteristics. Finally, correlations are presented with forced sinusoidal 
motion imposed about the pitch, yaw, and roll axes. For all force and moment coefficients, the reference area is 
Sref=0.77m 2 , c r e/=0.479m, and b re j=b/2=0.169m , MRP[x, y,z]=[0. 600, 0, 0]m. The pitch and yaw oscillations are 
imposed at the J RRP[x,y,z]=[0.855, 0, 0]m location, and the moments are transferred to the MRP. The support stings 
are not included in the force and moment integrations. As a reminder, the cases with undeflected control surfaces 
will be designated CSO, and the cases with deflected control surfaces LOB/LIB=-20/-20, ROB/RIB=+20/+20 deg as 
CS20. All moment coefficients are presented in the body axis system and are designated as Cmx, C M y, and C M z , for 
roll, pitch, and yaw, respectively. 

A. DLR-F19 static aerodynamics 

1. a- sweep 

A series of static time-accurate USM3D flow solutions were generated on the DLR-F19 for thirteen angles of 
attack between 0 and 30 degrees to explore the sensitivity of turbulence model on the medium grids G2.00 and 
G2.20 (see Tables 1 and 3). The SA and SST turbulence models were applied to both grids, and the k-8 model only 
to the G2.20 with control surface deflection CS20. 
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The CSO cases utilized a prescribed characteristic time-step of At* = 0.02 where physical At = At* - c re f /U^. 
Solutions were advanced in time for 2000 time-steps allowing the flow residual of the inner iterations to drop 2.5- 
orders of magnitude, which typically requires 5 to 6 inner iterations of the 2nd-order physical time-step scheme. 
This approach provided a total of 40 characteristic time steps for solution convergence, which was verified. The 
CS20 cases required a more aggressive strategy due to the additional complexity of asymmetric control surface 
deflection. Hence, a characteristic time-step of At* = 0.10 was prescribed for 2000 time steps, increasing the total 
number of characteristic time steps to 200. Some test solutions were run on the CSO case to confirm that the change 
in At* has negligible impact on the converged solution, which was expected due to the strong flow damping from 
Reynolds averaging. 

Comparisons of lift, drag, and pitching moment coefficient variations with angle of attack are presented in Fig. 8 
for both CSO (left) and CS20 (right) configurations. The flow physics behind the behavior of the experimental data 
[22] with increasing angle of attack, a, is well documented [11] and will not be repeated. The correlations of 
USM3D SA and SST turbulence model results for the undeflected CSO case is similar to that published during AVT- 
161 [24]. In general, both the SA and SST correlate well for C L and C D in the linear range up to a~15 deg. The C MY 
is under predicted in the linear range, which was typical of all RANS correlations from the AVT-161 task group [ 6- 
13]. A stronger sensitivity is observed at the higher post-stall angles. Note the large uncertainty bands on C MY for the 
experimental data above a>17 deg, which infer highly unsteady flow. The poor correlations of CFD in this range 
suggest that the inherent dissipation of the RANS formulation caused by excessive levels of viscosity damps any 
critical unsteady flow features that would influence this sensitive region of the angle-of-attack range. The SST 
model does yield a better correlation of post stall due to its known quality of producing lower levels of turbulent 
viscosity than SA for highly separated flows. 

The deflected CS20 correlations on the right of Fig. 8 reflect the general character of the CSO case. However, in 
addition to the SA and SST models, the KE (k-s) model has been added. The KE model offers no significant 
improvement to the correlations with experimental data. Furthermore, this model failed to run beyond cx=20 deg, 
thereby preventing any assessment in the post-stall region. 

2. /3- sweeps 

A series of static sideslip sweeps (Table 3) were generated on the DLR-F19 to explore the sensitivity of grids 
(Table 1) and turbulence models. The SA, SST, and KE turbulence models were applied to both CSO and CS20 
grids. 

A physical time step At=0.00167s was prescribed for the sideslip sweeps, which results in a characteristic time- 
step of At*= 0.165. A total of 25 inner iterations were performed between each time step. The sweeps were computed 
with one continuous run using the prescribe motion capability from FDCADRE. Restarting from a converged zero- 
sideslip initial solution at angle of attack, the configuration was moved in 1 degree increments of sideslip over 10 
time steps and then held on point for 290 time steps (At*~ 48) before moving to the next point. Convergence to 
steady state was confirmed within each increment. The advantages of this strategy are that it 1) requires a factor of 
3.5 less computer time than a comparable set of individual runs, and 2) mimics the wind tunnel test technique, which 
could potentially have path dependencies that affect the progression of leading-edge flow separation with increasing 
sideslip. 

The selection of the baseline grids for CSO and CS20 is determined from the results in Fig. 9, which convey the 
variations of static roll, pitch, and yaw moments, C M x, C M y, and C MZ , with /3 at the more challenging pitch angle 
0^=14.67 deg. The TN2545 experimental data [22] are plotted with symbols and uncertainty bounds. Note that the 
moments for the CSO case on the left are relatively symmetric, whereas those for CS20 with asymmetric control 
deflections are less so. Solutions were computed with USM3D/SA on grids Gl, G2, and G3 from Table 1. For the 
symmetric CSO cases, solutions were only computed in the positive /3-sweep direction (nose left). It is immediately 
obvious that Gl is not adequate for either CSO or CS20 cases. Grids G2 and G3 generally capture the slope of the 
moments (i.e. static stability), and to a lesser extent the inflection in Cmx and C M z at /3~±5 deg caused by a sudden 
change in flow state. It is also encouraging that grids G2 and G3 also produce the correct offset in negative rolling 
moment 0^-0.035 at /3=0 deg for CS20. Differences are observed in the offset and character of C M z for CS20, but 
most of the USM3D/SA solutions fall within the uncertainty range of the data for Cmx and C MZ . Considering the 
increasingly large cell count of grids G2 and G3, and their reasonably close correlation of results between these two 
grids in Fig. 9, the G2.00 (13 million cells) and G2.20 (18 million cells) grids will be selected for further assessment 
of turbulence models. 

Fig. 10 shows the effect of the SA, KE, and SST turbulence models on the static lateral/directional moment 
coefficients at the lower pitch angle 00=10.48 deg using the medium grids G2. Again, the CSO is on the left, and 
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CS20 on the right. The TN2545 experimental data, plotted with symbols and uncertainty bounds, exhibit a well 
behaved, almost linear quality. The USM3D solution results from the three turbulence models also correlate 
reasonably well with experiment. However, the SST model renders the best overall correlation. 

A similar correlation of turbulence models for the higher pitch angle <9o=14.67 deg in Fig. 11 yields a more 
surprising result. The SST model degenerates dramatically for both CSO and CS20 beyond (3 of ±2 degrees, whereas 
the SA and KE models retain reasonable correlations. The reason for this departure is evident in the corresponding 
surface C p contours and streamlines presented in Fig. 12 for fh+5 deg. At this sideslip, the rolling moment C M x in 
Fig. 11 for the SST solution on the CSO case shows a sudden departure in the positive direction away from the 
TN2545 data and the SA/KE solutions. Similarly the yawing moment C MZ for the CS20 solution with the SST model 
in Fig. 11 exhibits a similar sudden departure in the positive direction that extents well outside of the experimental 
uncertainty bounds. In Fig. 12, the SA and KE models produce a partially attached flow with accompanying low 
pressure along a majority of the right (windward) wing leading edge, and the beginnings of a leading-edge vortex 
near the left-wing tip. This is consistent with the production of a negative rolling moment Cmx . The SST model 
results in Fig. 12 reveal a well developed leading-edge vortex and accompanying low pressure on the left wing, and 
a more stalled flow with increasing pressures on the right wing, which is consistent with the production of a positive 
rolling moment Cmx . Numerous attempts were made to resolve the problem with the SST model, e.g. ran the finer 
G3 grids (Table 1), decreased time step At, increased number of inner iterations, decreased rate of incremental 
sideslip motion, initialized solution from free stream condition, etc., but nothing produced a significantly different 
result. In past experiences, the SST model has been a very reliable model of choice for many complex applications. 
However, it appears that the SACCON configuration at this higher angle of attack sideslip condition has exposed a 
shortcoming in the turbulence model. 

Overall in Fig. 11, the KE model yields the best correlations with TN2545 data. Unfortunately, KE solutions 
could not be obtained for /3 <- 6 deg. Hence, this model lacks the robustness for general use on the SACCON. For 
these reasons, the SA model is selected as the model of choice for SACCON S&C modeling, and will be utilized in 
the forced oscillation computations that follow. 

B. Forced oscillation (F-O) 

Sinusoidal oscillations about the pitch and yaw axes were imposed on the DLR-F19 SACCON at 
RRP[x,y,z]=[0.855, 0, 0]m through the FDCADRE trajectory input file. The roll oscillations for the DLR-F17 were 
prescribed by an internal motion generator within USM3D. The forced oscillation computations were initiated by 
solution restart from converged time-accurate static solutions at the nominal pitch angle 0 O (cfy), roll angle <P 0 = 0, 
and yaw angle W Q =0 deg. Although the dynamic force and moment coefficients are essentially converged to their 
periodic solution after the first quarter-oscillation cycle, the solutions are continued for two full sinusoidal cycles in 
this study. The second cycle is utilized for comparisons to eliminate any transients from the plots. The time- 
accurate RANS solution exhibits a deterministic quality by converging to a very stable solution, in contrast to the 
well known unsteadiness in the experimental data. The support sting moves with the wing in solid body motion, but 
is not included in the force and moment integrations. The moments are transferred from the RRP to the MRP. 

The USM3D convergence guidelines for dynamic RANS simulations are developed in Ref. [6]. In general, a 
well-converged solution is insured using a total of 36,000 iterations per oscillation cycle, regardless of the 
combination of outer time steps and inner iterations. This guideline is overly conservative for many conditions, but 
does insure convergence of the more difficult nonlinear cases using a single strategy. 

A physical time step of At=0.00167s was selected for the current computations to match that in the experimental 
data. For CSO cases with oscillation frequencies of 1.00Hz and 2.00Hz, this results in a At*= 0.174 and respective 
time-steps per cycle of 600 and 300 with corresponding inner iterations of 60 and 120 to maintain a total of 36,000 
total iterations per cycle. Similarly for CS20 cases with oscillation frequencies of 0.94Hz and 1.88Hz, the At*= 0.165 
with respective time-steps per cycle of 638 and 320 with corresponding inner iterations of 56 and 113. 

7. Pitch F-O ofDLR-F19 SACCON 

The cases for pitch oscillation of the DLR-F19 are listed in Table 4. Low and high frequency oscillations of 
A@=±4J deg about an initial pitch angle of <9o=10.3 deg are investigated for k P =2jif c re f/U^ of 0.06 and 0.12. These 
reduced frequencies correspond respectively to 1.00Hz and 2.00Hz for CSO, and 0.94Hz and 1.88Hz for CS20 cases. 

An important inconsistency in configuration was noted in Section II.B where the experimental TN2545 pitch 
oscillations [23] were imposed using a pitch link support with a 6° crank angle and an axis push rod. Hence, there 
was relative motion between the stationary post support sting and the oscillating SACCON configuration. For 
current computational results, the SACCON is supported by the 15° cranked yaw-link sting and pitch oscillations 
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were imposed as a solid body rotation of the entire grid, which contains no relative motions between the sting and 
wing. 

The pitch oscillation correlations are shown for the low frequency (£p=0.06) in Fig. 13, and high frequency 
(£p= 0.12) in Fig. 14. These figures portray the effect of grid density G1 and G2 using the USM3D/SA model on the 
aerodynamic response to the motion, which manifest as hysteresis loops for C N , C M y , and C M z that oscillate in pitch 
in the counter-clockwise direction between 5.6 < a< 15 deg. The effect of grid density is noticeable, but not large. 
The choice of turbulence model has a much greater impact as demonstrated in Ref. [24]. For computing dynamic 
stability derivatives, the shape of the hysteresis curves is the indicator of dynamic stability characteristics. In 
particular, the width of the hysteresis loop ACm=(C m + - C M ) at a 0 = 10.3 deg is proportional to the dynamic damping 
coefficient [2]. Here, C M + is the moment coefficient across a 0 with increasing angle, and C M - is the moment 
coefficient across a# with decreasing angle. (The vertical offsets observed between experiment and CFD are related 
to the static offset that often occurs from turbulence model or grid effects; these have no bearing on the dynamic 
stability derivative.) For this case the AC M y or AC M z at a 0 =@o = 10.3 deg are presented in Table 7. Because of the 
counter-clockwise direction of pitch and yaw hysteresis, the increment in moments are negative which indicate 
damping. In general, the pitch damping from CFD tends to be slightly higher than experiment, but is still reasonably 
close. The yaw damping due to pitch oscillation is essentially negligible. 


Table 7 - AC M y and AC M z coefficients for DLR-F19 undergoing pitch oscillation at O 0 =10.3° 





cso 



CS20 


kp 

AC m 

TN2545 

USM3D/SA 

TN2545 

USM3D/SA 




G1.00 

G2.00 


G1.20 

G2.20 

0.06 

AC my 
ACmz 

-0.007 

0.000 

-0.0097 

0.0000 

-0.0092 

0.0000 

-0.0062 

0.0002 

-0.0073 

-0.0004 

-0.0065 

-0.0004 

0.12 

AC my 
ACmz 

-0.013 

0.000 

-0.0169 

0.0000 

-0.0160 

0.0000 

-0.012 

0.0002 

-0.0143 

-0.0003 

-0.0123 

-0.0004 


2 . Yaw F-0 ofDLR-F19 SACCON 

The cases for yaw oscillation of the DLR-F19 are listed in Table 5. Low and high frequency oscillations of 
AW =± 5 deg about initial angles of O 0 =\ 0.1 and W Q =0 deg are investigated for kY=2jTfb re /Uoz of 0.10 and 0.19. These 
reduced frequencies correspond respectively to 1.00Hz and 2.00Hz for CS0, and 0.94Hz and 1.88Hz for CS20 cases. 
The SACCON configuration was supported by the 15° cranked yaw-link sting for both the wind tunnel [23] and 
CFD models with no relative motion between sting and wing during yaw oscillations. 

The yaw oscillation correlations are shown for the low frequency (£y=0.10) in Fig. 15, and high frequency 
(kr= 0.19) in Fig. 16. These figures portray the effect of grid density G1 and G2 using the USM3D/SA model on the 
aerodynamic response to the motion, which manifests as hysteresis loops for Cmx , C M y , and C M z that oscillate in the 
clockwise direction between -5 < W< +5 deg. As with the pitch oscillation, the effect of grid density for yaw F-0 is 
observed to be nominal. 

The width of the hysteresis loops AC M x and AC M z, where ACm=(C m + - C M ~) at Wg=0, are collected in Table 8. 
Here, C M + is the respective moment coefficient across W 0 with yaw angle increasing toward nose right, and C M - is the 
moment coefficient across W 0 with yaw angle decreasing toward nose left. The levels of ACmx and AC M z are very 
small due to the absence of any vertical surfaces. While very small, the ACmx does change sign from negative to 
positive between CS0 and CS20, indicating a tendency to lose roll damping. In general, the positive increments 
correspond to clockwise oscillation loops that represent propelling, and negative increments correspond to counter- 
clockwise loops that indicate damping. In general, the AC M from CFD is approaching 3 -times larger than experiment 
for the yaw F-O. While no specific reason has been discerned, this comparison is between two small numbers. 
Hence, in the absence of any vertical control surfaces and the small projected area of wing thickness these 
increments are highly sensitivity to the smallest perturbations. In contrast, the previous pitch and subsequent roll 
F-0 results experience a strong forcing function from the interactions of larger aerodynamic surface projections with 
motion. 
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Table 8 - AC mx and AC M z coefficients for DLR-F19 undergoing yaw oscillation at Q 0 =10.1° 





cso 



CS20 


ky 

ac m 

TN2545 

USM3D/SA 

TN2545 

USM3D/SA 




G1.00 

G2.00 


G1.20 

G2.20 

0.10 

AC mx 
ACmz 

0.000 

0.0003 

-0.0007 

0.0008 

-0.0007 

0.0009 

0.000 

0.0003 

0.0003 

0.0010 

0.0005 

0.0011 

0.19 

AC mx 
ACmz 

0.000 

0.000 

-0.0010 

0.0017 

-0.0010 

0.0018 

0.000 

0.0007 

0.0004 

0.0019 

0.0008 

0.0021 


3. Roll F-0 of DLR-F1 7 SACCON 

The cases for roll oscillation of the DLR-F17 with aft sting are listed in Table 6. Solutions were computed for 
seven frequencies of roll oscillation between A<P=±5 deg at angles of attack of a 0 = 10 and 20 deg. The reduced 
frequencies k R =2jtfb re f)/U 00 ranged between 0.063 and 0.264 (f= 0.24Hz to 1.00Hz). For sake of brevity, only three of 
these cases are presented in Figures 17 to 19 to document the roll F-0 hysteresis behavior and degree of correlation 
with experiment in a manner consistent with previous pitch and yaw results. Reference [25] utilized these solutions 
to explore the creation of reduce-order models from CFD data using system identification techniques. 

The plots of Cmx, C M z , and C Y vs. A<£> presented in Fig. 17, 18, and 19 correspond to k R = 0.063, 0.145, and 0.264 
(/= 0.24, 0.55, and 1.00Hz), respectively. The left column portrays the a 0 = 10 deg results, and the right column a 0 =20 
deg. For each frequency, dynamic damping is evidenced in roll and yaw from the counter-clockwise oscillation 
loops, and propelling (non-damping) for the side force as evidenced by clockwise loops. It should be mentioned that 
some vertical offsets present in the experimental data at zero roll angle, typically caused by small asymmetries in the 
model geometry and possible flow angularity, have been subtracted out [4] to facilitate more consistent comparisons 
of shape with the CFD predictions. The vertical offset is computed as the average of the two zero-axis crossing 
values. The T134 wind tunnel data in Fig. 17-19 also include variation bars denoting the standard deviation of flow 
unsteadiness. 

The effect of the SA/SST turbulence models on the coarsest grid Gla from Table 2 is shown in Fig. 17, 18, and 
19. In general for a 0 = 10 deg the CFD oscillation loops are in close agreement with the relatively steady 
experimental data and the effect of turbulence model is small, but less so at a 0 = 20 deg. Referring back to the static 
data in Fig. 8, a=10 deg corresponds to a more well behaved pre-stall condition, whereas a=20 deg falls within a 
very unsteady post-stall flow condition. Excellent correlations are exhibited Figures 17 to 19 at a= 10 deg for Cmx 
across the frequency range, correctly predicting the increase in roll damping with increasing roll frequency. While 
the trends with increasing frequency are reflected correctly for C MZ , and C Y at a 0 = 10 deg, there is a slight over 
prediction of the dynamic effect. At a 0 =2 0 deg the experimental data displays significantly more unsteadiness. As 
might be expected, the CFD solutions produce a more varied correlation with the data. An impact of turbulence 
model can also be observed at a#=20 deg with the SA model being in closer agreement with the mean value of the 
experimental data than the SST model. 

The effect of grids Gla/G2a/G3a from Table 2 using the SA turbulence model at a 0 =20 deg are investigated in 
Fig. 18 (right) for the mid- frequency casey=0.55Flz. While grid refinement does produce better correlations with the 
experimental data, the dominate effect comes from the SST turbulence model. An investigation of the impact of 
these CFD solution uncertainties on the creation of reduced order models constructed by system identification 
techniques is reported in Ref. [25]. There, two methods are presented, harmonic analysis and a two-step linear 
regression, for estimating mathematical models useful for stability and control analysis that can be applied to both 
CFD simulations and wind tunnel measurements. This approach offers a model structure that can accommodate 
unsteady aerodynamic behavior of the DLR-F 1 7 and characterize model parameter uncertainty. Regression analysis 
of the CFD/SST simulated data resulted in a steady-flow damping parameter in statistical agreement with the wind 
tunnel and CFD/SA cases but not for the unsteady term. For the CFD/SST case, the unsteady term was statistically 
different and with a much larger uncertainty bound than either wind tunnel or CFD/SA cases. Time history 
comparisons in Ref. [25] confirmed the reduced-order model agreement between the wind tunnel and CFD/SA 
model but the CFD/SST case showed the added uncertainty due to turbulence model selection. 
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V. Concluding Remarks 

The state of aerodynamic stability and control predictive capability of time-accurate RANS methodology has 
been investigated for a 5 3 -degree swept UCAV configuration with control surface deflections. The vehicle 
aerodynamics is dominated by the complex non-linear physics of round leading-edge vortex flow separation. 
Correlations with wind tunnel data are made for static longitudinal/lateral sweeps, and at varying frequencies of 
prescribed roll/pitch/yaw sinusoidal motion for the vehicle operating with and without control surfaces. The ultimate 
purpose for these investigations is to provide data and insight to researchers who are developing techniques for 
augmenting flight simulation models from low-speed CFD predictions of aircraft traversing non-linear regions of a 
flight envelope. 

These investigations included assessments of grid and turbulence model sensitivities with control surfaces 
deflected asymmetrically and undeflected. The effect of control surface deflection on the 3 -axis moment coefficients 
were modeled with comparable accuracy to that with undeflected surfaces for angle-of-attack and sideslip sweeps. 
For the static sideslip cases, the SST turbulence model exhibited premature leading edge flow separation and was 
deemed not suitable for further investigation. The SA turbulence model was selected as the most consistent approach 
for the full range of the test matrix. The roll/pitch/yaw forced oscillation computations produced dynamic stability 
characteristics consistent with that from the experimental data. 
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0.8406m 



Figure 1. SACCON planform and geometric parameters, controls surface designations, and location of 
surface pressure orifices (red dashed lines). 



Figure 2. DLR-F19 with deflected control surfaces mounted on 15° yaw-link in DNW-NWB tunnel (left). 
DLR-F17 mounted on aft-roll sting in NASA LaRC 14X22-Foot subsonic tunnel (right). 



right: G2.00, G2.20, G2a. 
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Figure 5. Volume grid through right-wing control surface of DLR-F19 CS20. Left to right: G1.20, G2.20, 
G3.30. 



Figure 6. Volume grid slices perpendicular to leading edge through right-wing and deflected control surface 
for DLR-F19 CS20. Left to right: G1.20, G2.20, G3.20. 



Figure 7. Close-up of right wing outboard flap gap for DLR-F19 CS20. Top: surface triangulation. Bottom: 
volume grid cut through control surface. Left to right: G1.20, G2.20, G3.20. 
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MY 







Figure 8. Effect of turbulence models on DLR-F19 for static lift, drag, and pitching moment coefficient. CSO 
(left) and CS20 (right). M^O.14, Re cref =1.40 million. NASA USM3D, Grid 2. 
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Figure 9. Effect of grid refinement on DLR-F19 static lateral/directional moment coefficients. CSO (left) and 
CS20 (right). M^O.14, Re crei =1.40 million, 00=14.67 deg. NASA USM3D/SA. 
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Figure 10. Effect of turbulence model on DLR-F19 static lateral/directional moment coefficients. CSO (left) 
and CS20 (right). M^O.14, Re crei =1.40 million, 00=10.48 deg. NASA USM3D, Grid 2. 


18 

American Institute of Aeronautics and Astronautics 



Figure 11. Effect of turbulence model on DLR-F19 static lateral/directional moment coefficients. CSO (left) 
and CS20 (right). M^O.14, Re cref =1.40 million, 00=14.67 deg. NASA USM3D, Grid 2. 
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Figure 12. Effect of turbulence model on DLR-F19 surface C p contours and streamlines at /3=+5 deg. 
CSO (left) and CS20 (right). IYl y =0.14, Re cref =1.40 million, ©,,=14.67 deg. NASA USM3D, Grid 2. 
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Figure 13. Effect of grid density on DLR-F19 longitudinal and directional coefficients under low frequency 
forced PITCH oscillation (k P =0.06), 00=10.3 deg, A0=±4J deg. CSO (left) and CS20 (right). M^O.14, 
Re crei =1.5 million. NASA USM3D/SA, Grids 1 and 2. 
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Figure 14. Effect of grid density on DLR-F19 longitudinal and directional coefficients under high frequency 
forced PITCH oscillation (k P =0.12), 00=10.3 deg, ^0=±4.7 deg. CS0 (left) and CS20 (right). Moo=0.14, 
Re cref =1.5 million. NASA USM3D/SA, Grids 1 and 2. 
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Figure 15. Effect of grid density on DLR-F19 lateral/directional moment coefficients under low frequency 
forced YAW oscillation (k Y =0.10), 00=10. 1 deg, A 1 F=±5.0 deg. CSO (left) and CS20 (right). M^O.14, 
Re cre i=1.5 million. NASA USM3D/SA, Grids 1 and 2. 
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Figure 16. Effect of grid density on DLR-F19 lateral/directional moment coefficients under high frequency 
forced YAW oscillation (k Y =0.19), @0=10.1 deg, A l P=±5.0 deg. CSO (left) and CS20 (right). M^O.14, 

Re cre i=1.5 million. NASA USM3D/SA, Grids 1 and 2. 
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Figure 17. Effect of turbulence model on DLR-F17 lateral/directional coefficients under forced ROLL 
oscillation f=0.24Hz, k R =0.063, AO=±5 deg. a=10 deg (left) and cx=20 deg (right). M^O.14, Re cre f=1.40 
million. NASA USM3D/SA/SST, Grid la. 
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Figure 18. Effect of grid and turbulence model on DLR-F17 lateral/directional coefficients under forced 
ROLL oscillation f=0.55Hz, k R =0.145, AO=±5 deg. cx=10 deg (left) and a=20 deg (right). Moo=0.14, Re cre f=1.40 
million. NASA USM3D/SA/SST, Grids la and 2a. 
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Figure 19. Effect of turbulence model on DLR-F17 lateral/directional coefficients under forced ROLL 
oscillation f=1.00Hz, k R =0.264, AO=±5 deg. a=10 deg (left) and a=20 deg (right). Moo=0.14, Re crel =1.40 
million. NASA USM3D/SA/SST, Grid la. 
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